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In the framework of an extended Nambu-Jona-Lasinio model we are studying pion condensation 
in quark matter with an asymmetric isospin composition in a gravitational field of the static Einstein 
universe at finite temperature and chemical potential. This particular choice of the gravitational 
field configuration enables us to investigate phase transitions of the system with exact consideration 
of the role of this field in the formation of quark and pion condensates and to point out its influence 
on the phase portraits. We demonstrate the effect of oscillations of the thermodynamic quantities 
as functions of the curvature and also refer to a certain similarity between the behavior of these 
quantities as functions of curvature and finite temperature. Finally, the role of quantum fluctuations 
for spontaneous symmetry breaking in the case of a finite volume of the universe is shortly discussed. 
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I ■ I. INTRODUCTION 

Low-energy nonperturbative effects in quantum chromodynamics (QCD), especially at nonzero densities, can only 
be studied by approximate methods within the framework of various effective models. It is well known that light meson 
physics is described by four-fermion models, such as the Nambu-Jona-Lasinio (NJL) model, which was successfully 
used to deal with dynamical chiral symmetry breaking (D^SB) both in the vacuum and in hot/dense baryonic matter 

i^T (see, e.g. [l|, [2[, as well as the reviews 0, El and references therein). Recently, much attention has been paid to the 

effects of diquark condensation and color superconductivity (CSC). The first studies of the gap equations and the 

f-^. . Ginzburg-Landau free energy for a system of relativistic fermions led to the conclusion that superconductive and color 
superconductive states may arise in baryonic matter [f|, Q (see also the recent papers p], |8|). Another interesting 

_il , phenomenon, the condensation of charged pions, which may appear in dense hadronic matter due to an asymmetry 
of its isospin composition, has been investigated in the framework of QCD-like effective models, including the NJL 
model, as well [|, M EJ 111 E M El ■ 

Note that all these phenomena might be inherent to physics of compact stars, where rather strong magnetic as well 
as gravitational fields are present. Therefore, investigations of the influence of an external (chromo- ) magnetic field on 
the properties of the D%SB phase transition [lfl [ljl. color superconductivity [l8| and pion condensation [19| effects are 
quite appropriate. In particular, it was shown in [16l.ll7l.ll8l.ll9l] that external fields significantly change the properties 
of the chiral and CSC phase transitions. In several papers, in the framework of the NJL model, the influence of a 
gravitational field on the D^SB due to the creation of a finite quark condensate (qq) has been investigated at zero 
values of temperature and chemical potential [2(J, |2l], [22|, [23| • The study of the combined influence of curvature and 
temperature has been performed in [241 ]. Recently, the dynamical chiral symmetry breaking and its restoration for 
a uniformly accelerated observer due to the thermalization effect of acceleration was studied in [251 ] at zero chemical 
potential. Further investigations of the influence of the Unruh temperature on the phase transitions in dense quark 
matter with a finite chemical potential, and especially on the restoration of the broken color symmetry in CSC were 
made in [2a |. 

One of the widely used methods of accounting for gravitation is based on the expansion of the fermion propagator 
in powers of small curvature R [27|, (28|. For instance, in [29|], the three-dimensional Gross-Neveu model in a spacetime 
with a weakly curved two-dimensional surface was investigated, using an effective potential at finite curvature and 
nonzero chemical potential. In paper [30J, this weak curvature expansion was used in considering the D^SB at non- 
vanishing temperature and chemical potential. It should, however, be mentioned that near the phase transition point, 
one cannot consider the critical curvature R c to be small and therefore the weak curvature expansion method can not 
be applied. Hence, in the region near the critical regime different nonperturbative methods or exact solutions with 
finite values of R should be used. This kind of solution with consideration for the chemical potential and temperature 
in the gravitational background of a static Einstein universe has been considered in [3l( . There it was demonstrated 



that chiral symmetry is restored at large values of the space curvature. Analogous studies of diquark condensation 
and the related color symmetry breaking under the influence of a gravitational field have been performed recently in 
the model of a static Einstein Universe in [33] ■ Recall that this model is widely discussed in literature either as a 
solution of the Einstein equations with a given cosmological constant and a nonvanishing energy-momentum tensor 
of an ideal fluid as a source, or as an initial state in inflationary cosmology with a scalar field, and the cosmological 
constant as its vacuum energy (see, for instance, [33j] ). Moreover, the Einstein universe and other suitably generalized 
compact curved spacetimes were extensively employed for studying the phenomenon of Bose-Einstein condensation 
(see, e.g. [34J and references therein). As is well known, one of the possible models for the expanding universe is the 
closed Friedmann model. Since the formation of quark condensates is expected to take place considerably faster than 
the expansion rate of the Universe, its radius can be considered in our calculations as constant. In this sense, the 
chosen model of a static Einstein universe can be considered as a simple cosmological model for a space with positive 
curvature. On the other hand, the same form of the metric does also hold for the interior of a (collapsing) spherically 
symmetric star, thereby admitting also a non-cosmological interpretation of the chosen model. 

As was mentioned above, in quark matter there may arise another remarkable phenomenon, pion condensation. 
Up to now little information about the properties of pion condensation in the presence of external fields has been 
obtained. In particular, there arises the interesting question, what influence gravitational fields can produce on pion 
condensation. The main aim of the present paper is just to give a detailed investigation of this issue. 

In order to be able to consider the effects of gravitation on pion condensation in a rigorous nonperturbative way, 
we will investigate here the phase structure of isotopically asymmetric quark matter in the framework of an extended 
NJL model with dynamical breaking of chiral and flavor symmetries in the static Einstein universe. In particular, our 
calculations are performed in the most simple case, restricting ourselves, for simplicity, to the flavor 517(2) group. 
By using the mean field approximation, we will then derive an analytical expression for the thermodynamic potential 
of the NJL model that will enable us to calculate the chiral and pion condensates of quarks moving in a space with 
constant positive curvature. On this basis, we study, in particular, the phase portraits of the model. Our paper 
is organized as follows. Section II, III contain the Nambu-Jona-Lasinio model in curved spacetime and the general 
formalism for the derivation of the thermodynamical potential. Numerical calculations and conclusions are given in 
sections IV, V. Finally, some estimates of the role of quantum fluctuations in the case of a closed universe are given 
in the appendix. 

II. NAMBU - JONA-LASINIO MODEL IN CURVED SPACETIME 

Suppose that dense, isotopically asymmetric quark matter (in this case the densities of u and d quarks are different) 
in curved spacetime is described by an extended NJL model with the following action: 



and the Lagrangian 



C-i 



(1) 



i7"V„ - m + /i7° + 5fiT 3 j° q + G {qqf + (qij 5 Tq) 2 , (2) 



where the quark field q{x) = qi a {x) is a flavor doublet (i = 1,2 or i — u,d) and color triplet (a = 1,2,3) as well as 
a four-component Dirac spinor (the summation in ([2]) over flavor, color and spinor indices is implied); t% (i — 1,2,3) 
are Pauli matrices. In 4-dimensional curved spacetime with signature (+, — , — , — ), the line element is written as 

ds 1 =T ]hl e\ i eldx^dx v . 

The gamma- matrices 7^, the metric g^ v and the vierbein e a , as well as the definitions of the spinor covariant derivative 
V v and spin connection u>^ h are given by the following relations [28j : 

{l l i{x)^ v {x)} = 2g tiV (x), {7a,7g} = 2r/ aS , r? aS = diag(l, -1,-1,-1), 

9^9 vp = SP, g^(x) = e^(x)e v& (x), f M (x) = efa)<ft. (3) 

V„=d fl + T M , r M = - wf a a , a hh = -[7a, 7g], 

w ju — 2 e ° e p ^P^pv- - CpXf, . - C^p], Cxpfj, — eld[ p e^a- (4) 



Here, the index a refers to the flat tangent space defined by the vierbein at spacetime point x, and the 7 a (d = 0, 1, 2, 3) 
are the usual Dirac gamma-matrices of Minkowski spacetime. Moreover 75, is defined, as usual (see, e.g., [28|]), i.e. to 
be the same as in flat spacetime and thus independent of spacetime variables. 

In order to describe the quark composition of matter we introduced fi = fis/3 in ^ as the quark number chemical 
potential. Since the generator 73 of the third component of isospin is equal to T3/2, the quantity 5fj, in ([2]) is half the 
isospin chemical potential, 5fi = fxi/2. If the bare quark mass m is equal to zero, then at Sfi — 0, apart from the 
trivial color SU(3) symmetry, the Lagrangian @ is invariant under the chiral group SUl(2) x SUr(2) transformations. 
However, at S^j, 7^ this symmetry is reduced to the subgroup C/j 3 l(1) x Ui 3 r(1) (here and above the subscripts L, R 
imply that the corresponding group acts only on the left, right handed spinors, respectively). It is convenient to 
present this symmetry as Z7j 3 (l) x Uai 3 {1), where Uj 3 (l) and Uai 3 {\) are the isospin and axial isospin subgroup of 
the chiral SUl(2) x SUr(2) group. Quarks are tranformed under these subgroups in the following way 



tfr,(l) 



exp(iar 3 )g, 



Uai 3 {1) 



exp(ia 7 r 3 )q, 



(5) 



where a, a' are independent group parameters. At nonzero bare quark mass, m ^ 0, and nonzero isotopic chemical 
potential, i.e. Sfi 7^ 0, the Lagrangian ^ is still invariant under the isospin subgroup t// 3 (l), but invariance with 
respect to Uai 3 (1) does no longer hold in this case. 

The linearized version of Lagrangian @, that contains collective bosonic fields cr(x) and iTk(x) (k — 1, 2, 3), has the 
following form 



C 



1 
4G 



J il v ^u + M7° + Sfir 3 j a - a - m — i7 5 7r fc r fc 

From the Lagrangian ([6]), one can find the equations of motion for the bosonic fields, 

a(x) = -2G(qq); n k (x) = -2G(gi 7 5 r fe g). 



TTfcTTfc 



(6) 



(7) 



It is clear from these relations that under isospin Uj 3 (l) and axial isospin [7aa 3 (1) transformations (J5J the bosonic 
fields ([7|) are changed in the following way: 



Uj 3 (l) : a —> a; n 3 — > 7r 3 ; 7Ti — ► cos(2a)7r! + sin(2a)7r 2 ; 7r 2 — > cos(2a)7r 2 — sin(2a)7Ti, 
7 7 J 4/ 3 (1) : m — >• 7Ti; 7r 2 — > 7r 2 ; cr — > cos(2a')cr + sin(2a')7r 3 ; 7r 3 — » cos(2a')7r 3 — sin(2a')a. 



(8) 



In the fermion one-loop (mean field) approximation, the effective action for the boson fields is expressed through 
the path integral over quark fields: 



exp(i<S c ff (<t, 7Tfe)) — N' I [dq] [dq] exp ( i / d 4 XyJ—g C 



where 



<Seff(cr,7T fc ) = - / d A x^J~^g 



a 2 +irl 



AG 



- S c ft, 



(9) 



N' is a normalization constant. The quark contribution to the effective action, i.e. the term S g in ((9|), is given by: 

exp(j,5cff) = N' [dq][dq] exp(i / d 4 xy/^g qVq) = dot V. (10) 

In (jTUJ) , we have used the following notation 



V = vfV v + M7° + <5/-iT37° - o - m 



IT^kTk- 



(11) 



Clearly, T> is an operator in the coordinate, spinor and flavor spaces. Apart from this, it is proportional to the unit 
operator l c in the color space. Thus, it can be presented in the flavor<g)color space in the following matrix form: 



V 



A , B 

B , A 



1, 



d® a r . 



(12) 



./■ 



where operators A, A, B, B act in the coordinate and spinor spaces only, and 

A = ij u V„ + (/i + <5/i)7 — a — m — ij J ir 3 , B = 



A = i'fVu + (/i - 5/i)7° 



cr 
cr 



Z7 5 (7Ti 



Z7T 2 ), 



m-M7 5 7r 3 , 73 = i7 5 (7Ti + Z7r 2 ). 



(13) 



Due to the trivial color structure, it follows from (fT2l) that 



det V = (det D) 3 . (14) 

Now, suppose that ttt, — n2 — and cr, tt\ are quantities independent of coordinates . In what follows, we shall denote 
the pion condensate ~k\ by A for convenience. Then, using the general formula 

det ( g ' B A J = dct[-BB + BAB' 1 A] = det[AA - ABA~ l B], (15) 

we obtain 

r rcr 2 + A 2 l f 

S cS {a, A) = - / d 4 x^/=g~ J - 3ilndetL> = -Q(a, A; /x, SfJ,) / d 4 x^/=g~, (16) 

where we have introduced the thermodynamic potential fi(cr, A;//, b\x) of the system at zero temperature and 

detD = detJA 2 + \-i-fX! v - + 6p)j° -a- m][i^V u + (fj, - 6fj,)-y° -a- m]\. (17) 

III. THERMODYNAMIC POTENTIAL 

A. General formalism 

The line element in the static Einstein universe is defined by the following relation: 

ds 2 = g fn ,(x)dx f "dx" = dt 2 - a 2 (d X 2 + sin 2 X (d9 2 + sin 2 8dtp 2 )), (18) 

where a is the radius of the Einstein universe (this quantity is related to the scalar curvature, R = 6/a 2 ); — oo < t < oo, 
0<x< 7r i0<^< 7r i0<ip< 2ir. Clearly, 7°(a;) in this case anticommutes with all j k (x) and commutes with all V„, 
where v = 0, 1, 2, 3 and k = 1, 2, 3. Moreover, Vo = <9o- So, starting from (fTT)) . we have (^ = 0, 1, 2, 3; k, I = 1, 2, 3): 

det£> = dctJA 2 + ^V,,) 2 - i(/i + ^hV V y - i(/i - (5^)7 ,y V l/ 7 - fi 2 + Sfi 2 + 26/J,(a + m) 7 ° + (cr + mf\ 

= detJA 2 + d 2 - 2i[id a - /i 2 + 7 fe V fe V V; + (cr + m) 2 - 25(i[i"? -f k V k - (cr + 7/1)7°] + <5/i 2 }. (19) 

Let P = i9 , U = -i7°7 fc V fe + (ct + 771)7°. It is evident that TCH = 7 fc V fe 7'V; + (cr + tti) 2 . Hence, we have from (fT9|) 

det£> = detO, where = A 2 - (P Q + fi) 2 + (H - 5[i) 2 . (20) 

Evidently, H is an operator in the Hilbert space of functions, depending on the space coordinates x. As it is 
well-known (see, e.g. [3a, [36|), each of the eigenvalues ±Ei of this operator is c^-fold degenerate, 

Ei = v /cj 2 + (m + cr) 2 , w , = i(/ + |), d, = 2(i + l)(J + 2), J = 0,1,2.... (21) 

Thus, one can write 

Hipi a n(x)=r]Eiipi ari (x), / d^x^/^g^ia^xjipva^^x) = 6w 6 aa > 6^ , (22) 



where the eigenfunctions ipi ar] {x) (a = l,...,di;rj = ±1) of the operator 7i are also known (see, e.g., [35l l36j|). and 
9 = det g M „ = —det gij(x) (see (fl"8|l ). Now, let us choose in the Hilbert space of functions a basis of the form 
^lar/po (t, x) = e~ lpot ^i a r)(x), where — oo < po < oo. Since each element ^i ariPo (t,x) of this basis is an eigenfunction 



1 To justify this assumption, let us consider for simplicity the case m = 0. If the fields a, ir^ do not depend on coordinates, then the 
effective action is a function of the invariants c 2 + 7r| and 7r 2 + 7r| only (this fact is due to the symmetry of the model with respect to 
the transformations J5}). Therefore, without loss of generality one can put 773=712 = 0. 



both of Pq and TC, one can easily conclude that the operator O from ([20]) is diagonal in this basis, i.e. each ^i aV p {t, x) 
is an eigenfunction of O. The corresponding eigenvalues £-i r)Vo of O have the following form: 

£ lvPo = A 2 - (po + ^) 2 + (vEi + Sfi) 2 . (23) 

It is clear from (|23j) that eigenvalues £z,jp of the operator O do not depend on the quantum number a = 1, ...,d;, 
being c^-fold degenerate. Taking into account in (fT6|) the relations (|20|) and DetO = exp(TrlnO), as well as the results 
of Appendix [Al where the quantity TrlnO is calculated (see (JA11I) ). one finds 

S cB (a, A) + / d^x^-g l ^tf ] = -3zTrlnO = -3zT £ / ^°d« In [A 2 - (p + ^) 2 + (^ + <^) 2 ] • (24) 
^ L J lr] J 7r 

Here J d 4 Xy/—g = TV, where T = J dt stands for an infinite time interval and V = J d 3 Xy/—g = 2ir 2 a 3 is the space 
volume of the Einstein universe (the last relations are due to the fact that g = — dct gij(x) depends only on x). Now, 
using the definition (|1€>[) of the thermodynamic potential (TDP) and summing in (|24[) over rj = ±1, we have for the 
zero temperature case 

fi( CT) A ;/ ^) = ^^ + |f; /'^{ln[A 2 -( P0 + M ) 2 + (£; ( -^) 2 ]+ln[A 2 -b + / ,) 2 + (^ + ^) 2 ]}. (25) 



i=0 



To find the TDP fl(a. A; fi, Sfi, T) in the case of nonzero temperature T, one should use the imaginary time technique, 
where, after summation over Matsubara frequencies (see, e.g., |32|.|37|). the following expression can be found 

2 I A 2 Q °° 071 °° 



4G V<^ L ' ' J V 

2=0 z=o 



3T 
V~ 



^d / {ln[l + e-^" )+ ' i )] +ln[l + e -^ £ ^ ) -^]} (26) 



with E) = \J [Ei ± Sfi) 2 + A 2 and /3 = 1/T. It is clear that Q(a,A;/i,Sfi,T) is an even function with respect to 
each of the transformations fi — ► — jjl or Sfi — ► —Sfi. Thus, one can deal only with non-negative values of the chemical 
potentials, fi > 0, Sfi > 0. 

From this moment on, we will consider only the case of nonzero isospin chemical potential Sfi 7^ 0, whereas the 
baryon chemical potential is set equal to zero, fi = 0, since its presence is not of principle importance for us. So, at 
fi = two particular cases can be investigated on the basis of the TDP (f26|) . 

First, let us choose T = and fi = 0, but Sfi 7^ 0. Then we obtain from (|26|) the expression: 

flfo A; M = 0,<5m, T = 0) = fi(<r,A;fy) = ^^_! - i^d, {£<+> +^ ( " ) }. (27) 

2=0 

Secondly, at T 7^ 0, /i = 0, Sfi 7^ we obtain: 

fifo A;m = 0,<5/x,T) = fi(er, A;8fi,T) = ^— - tj 2>{25, {+) +4~ ] } 

1=0 

-^-^^{ln[l + e -^ <+) ]+ln[l + e-^^]}. (28) 

2=0 

Next, let us consider the limit of zero curvature or infinitely large radius of the universe. It is clear that the metric 
(|18p never coincides with that of the flat Minkowsky spacetime because these two spacetimes have different topologies. 
However, in the limit a — > 00 and R — > one can obtain from (|28|l the usual expression for the TDP in flat spacetime 
(see, e.g., [lj]) by the following substitution: 



k, toi^k, d l = 2(l + l)(l + 2)^2k 2 a 2 , ^ -> / dl = a j < 



In this case, the TDP looks as follows: 



<t 2 + a 2 „ r d 3 k 



, f d 6 k ; 



Jl( t r J A;*/i ) T) = ^^=--6 / 7^ ] 4 +) + 4 _) + 2Tln[l + e"^ +) ] + 2Tln[l + e"^"']} , (29) 



where E { k ±] = y^ ± 5fi) 2 + A 2 . 

It should be noted that one more particular case, when T 7^ 0, \i 7^ 0, but 8fx = 0, can easily be reduced to the 
investigation of the one-flavored NJL model at T 7^ 0, (X 7^ in the Einstein universe [31(, and hence we shall not 
consider it here. 

B. Regularization 

First of all, in order to normalize the TDP, we should subtract a corresponding constant from it, such that f2(cr — 
0, A = 0) = 0. The thermodynamic potential, normalized in this way, is still divergent at lar ge I , and hence, we 
should introduce a (soft) cutoff in the summation over I by means of the multiplier e~ Ul ' A [3ll |32|. where A is the 
cutoff parameter 2 . 

For convenience, we shall multiply all dimensional quantities that enter the thermodynamic potential by the corre- 
sponding power of A to make them dimensionless, i.e., Vl/ A 4 , er/A, A/A, A 2 G, A 3 V, i?/A 2 , T/A, ///A, Sfi/A, oJi/A, 
and denote them using the same letters as before: 51, a, A, G, V, R, T, /j,, 6fi, wj. Then the regularized thermody- 
namic potential is written as 

n^(a,A;S„,T) = *±* if] e -"d,{i5< + > + E^} 



AG V , 

1=0 



-^Y,e- bJl di{\n[l + e-^ +) } + In [l + e"^"'] }. (30) 

1=0 

In the following Section, we shall perform a numerical calculation of the points of the global minimum of the finite 
regularized thermodynamic potential Jl rcg (cr, A) — i! rcg (0, 0) (they should of course coincide with the minima of the 
potential f2 rog (er, A)), and with the use of them, consider quark matter phase transitions in the gravitational field of 
the Einstein universe. 

IV. NUMERICAL CALCULATIONS 



In this section, on the basis of the thermodynamic potentials (127l) - (|28[) . we will study numerically phase transitions 
in quark matter and consider only the case of nonzero isospin chemical potential Sfi ^ 0, whereas the baryon chemical 
potential is set equal to zero, /1 = 0. In order to obtain the values of condensates, one should find the global minimum 
point (GMP) of the thermodynamic potentials over the variables a, A from the corresponding gap equations 

dn re s(a,A) _ dn rc s(cr, A) _ 

Formally, there are four different expressions for the GMP: i) (a = 0, A = 0), ii) (a 7^ 0, A = 0), iii) (a = 0, A 7^ 0), 
iv) (cr 7^ 0, A 7^ 0). The first two GMPs correspond to the isotopically invariant phases of the model, whereas the 
GMPs of the form iii) and iv) correspond to the phases, in which the ground state is no more £// 3 (l) invariant. In 
these phases the pion condensation phenomenon occurs. For simplicity, we take for the numerical calculations of the 
GMP the value of the coupling constant G = 1 (in our dimensionless choice of parameters). 



2 In flat spacetime, the cutoff constant A can be specified according to the experimental results. However, in the case of a curved spacetime, 
in order to fix the cutoff A, we need new theoretical/experimental inputs for chiral QCD in (strong) gravitational background fields, 
concerning, for instance, an effective gluon mass, known values of the quark condensate or even experimentally measured characteristics 
of a pion. Due to the lack of experimental knowledge, our aim is here to perform only a qualitative study of gravitational effects on the 
quark and pion condensates by investigating the respective phase portraits of the NJL model. For this reason it is convenient to scale 
the thermodynamic potential and all relevant quantities like condensates, curvature, chemical potential and temperature by the cutoff 
A. 



A. Zero temperature 

Let us first consider phase transitions at zero temperature, T — and choose the current quark mass to be equal 
to zero, m = 0. The thermodynamic potential in this case is described by formula (127|) . 

The detailed investigation of the GMP properties vs external parameters R and Sfi results in the phase portrait 
shown in Fig. Q] For the points of symmetric phase 1, the GMP is at a = and A = 0. In the phase 3, the minimum 
is at a = and A ^ 0, and this indicates that the isospin Uj 3 (l) symmetry of the model is dynamically broken in 
this phase. We note that at the same time the chiral Uai 3 (1) symmetry remains unbroken in this phase, since in the 
GMP we have a = 0, as it should be in the case of zero current quark mass, when A^O. 



<5ju 




R 



10 20 30 40 



FIG. 1: The phase portrait at zero temperature T = and m = 0. Number 1 denotes the symmetric phase and 3 denotes the 
isospin symmetry breaking phase (the phase with the pion condensate A/0). 
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FIG. 2: The behaviour of the pion condensate A. Left picture: T — 0, 5/i — 4.5. Right picture: R = 15, T = 0. 



As one can see from Fig. [TJ the critical curve, which separates the phases 1 and 3, has an oscillating character. 
This phenomenon is explained by the discreteness of the fermion energy levels (|21[) in compact space. Moreover, 
as it is clearly seen from the curves in [5J the pion condensate A vs R also oscillates in the phase 3. This effect 
resembles the van Alphen-de Haas oscillations of different phy sical quantities in the magnetic field H, where fermion 
levels are also discrete (the Landau levels) [38| (see also [3!|, where a similar influence of a magnetic field on the 
oscillation behavior of the Compton scattering and photoproduction cross-sections was demonstrated). Indeed, the 
corresponding magnetic oscillations of the critical curve in the \i-H phase portrait of dense cold quark matter with 
four-fermion interactions were found in papers [161 ]. There, the existence of the standard van Alphen-de Haas magnetic 
oscillations of some thermodynamical quantities, including magnetization, pressure and particle density of cold dense 
quark matter was also demonstrated. 
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The behavior of the pion condensate A as a function of R at fixed 5/j, and as a function of 5(i at fixed R is shown 
in Fig. [2] (left and right pictures respectively). 

The phase portrait at finite current quark mass, m 7^ 0, is depicted in Fig. [3l In phase 2 the chiral symmetry is 
now broken due to a finite value of the current quark mass, and the global minimum of TDP is at a 7^ and A = 0. 
In the mixed phase 4 both condensates are nonzero, i.e. a 7^ and A/0. 



S/.1 



4 

2 






10 



20 



30 



40 



R 



FIG. 3: The phase portrait at zero temperature, T = 0, and m = 0.01. Number 2 denotes the chiral symmetry breaking phase 
with u/0,A = 0, and 4 denotes the mixed phase with a 7^ and A 7^ 0. 

The behavior of the chiral condensate a in the case of finite quark mass, m 7^ 0, is shown in Fig. [4] as a function of 
R at 5(i = 4.5 (left picture) and as a function of 6(1 at R = 15 (right picture). 




LT 




FIG. 4: Condensate a at 5(i = 4.5 (left picture) and R — 15 (right picture), m = 0.01 

One can see oscillations of a on both pictures, although they are rather strong in the left picture for the dependence 
on R, while they are weakly seen in the high 5/j, tail of the curve in the right picture. One should note that the right 
picture resembles, except for these oscillations, the corresponding curve in [11] (Fig. 1) for the flat case. 



B. Finite temperature 



Using formula (|30[) for the thermodynamic potential, we can also study the influence of finite temperature on phase 
transitions. The phase portrait at T = 0.1 and zero current quark mass, m — 0, is shown in Fig. [5] in terms of R— 5fi. 
It is seen from this figure that growing temperature leads to a smoothing of oscillations of the phase curve. 

For comparision, in Fig. O the phase portraits in the T — 5(i and R — 5(i planes (the latter now at another value of 
temperature T = 0.4) are depicted. First of all, it is clear from Fig. [6] that the isospin symmetry is restored due to 
the vanishing of the pion condensate both at high temperature and high curvature. The similarity between these two 
plots leads to the conclusion that curvature and temperature effects play a similar role in the restoration of isospin 
symmetry. 




R 



FIG. 5: The phase portrait at T — 0.1 and m — 0. Number 1 denotes the symmetric phase and 3 denotes the isospin symmetry 
breaking phase (the phase with the pion condensate A ^ 0). 
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FIG. 6: The phase portraits at R = 4 (left picture) and at T — 0.4 (right picture), m — 0. 

The phase portrait in Fig. [7] for R = corresponds to the case of flat Minkowski spacetime. It looks very similar 
to the one obtained, for instance, in |ll( (see upper panel of their corresponding Fig. 12). Let us compare the phase 
portraits at zero (Fig. [7]) and finite curvature (Fig.[|5J left picture). In the first case the pion condensation appears at 
arbitrary small nonzero values of the isospin chemical potential, while in the second case the isospin symmetry becomes 
dynamically broken only at some finite value of the chemical potential <5/i. This phenomenon may be explained by 
the existence of a gap in the quark spectrum (|21[) . which is proportional to the inverse radius of the Einstein universe. 
In this case the effect of curvature is similar to the effect of a nonzero current quark mass (for comparison see lower 
panel of Fig. 12 in [llj and our Fig. [B]). 

V. CONCLUSIONS 



In the framework of an extended Nambu-Jona-Lasinio model, we have studied the influence of a gravitational field on 
pion condensation in isotopically asymmetric quark matter at finite temperature and isospin chemical potential. As a 
particular model of a gravitational field configuration we have taken the static Einstein universe. This particular choice 
enabled us to investigate phase transitions of the system with an exact consideration of the role of the gravitational 
field in the formation of the quark and pion condensates and thus to demonstrate its influence on the phase portraits. 
In particular, we have found that thermodynamic quantities such as quark and pion condensates as well as the 
corresponding phase boundaries (critical curves) oscillate as functions of curvature. This oscillating behavior is 
smoothed out with growing temperature. There exists also an interesting similarity between the behavior of phase 
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12 3 4 5 

FIG. 7: The phase portrait at R — and m — 0. 

portraits when considered as functions of curvature and as functions of finite temperature (compare Figs. 6). Moreover, 
we have shown that for massless quarks and for some values of R the isospin symmetry becomes in curved spacetime 
dynamically broken only at some finite value of the chemical potential Sfj, (look at Fig. 6 (left picture) for rather 
small values of T). This is contrary to the case of flat spacetime, where the pion condensation appears at arbitrary 
small nonzero values of the isospin chemical potential (see Fig. 7) . This effect resembles the pion condensation but 
at nonzero bare quark mass and may be explained by the presence of a gap in the energy spectrum of quarks in the 
gravitational field. 

Finally, let us add some remarks on the possible role of quantum fluctuations of the collective fields and finite 
size effects in our approach. In fact, since the volume of the space region modelled by the closed Einstein universe 
adopted in this paper is limited, finite size effects could eventually change the character of phase transitions (see e.g. 
discussions in [40, El]]). Eventually, this might even lead to particular situations, where no phase transition can occur. 
However, as it is clear from physical considerations, the finite size in itself may not practically forbid the dynamical 
symmetry breaking, if the characteristic length of the region of space occupied by the system is much greater than 
the Compton wavelength of the excitation responsible for tunneling and restoration of symmetry (see, e.g., |4fj|). It 
should be noted that our results, obtained with the use of the mean field approximation in the framework of the NJL 
model are consistent with corresponding estimates of the role of quantum fluctuations (see Appendix B). 

In conclusion, we emphasize that the results of this paper are evidently only of a qualitative nature. They do 
not allow us to find an exact value of the critical radius, and hence, further studies with more realistic models of 
gravitational fields should be undertaken. 
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APPENDIX A: THE TRACE OF THE OPERATOR 6 (f20|) 

In Section IIIII we have introduced two operators, Po and H, acting in the Hilbert space H of all quadratically 
integrable functions f(t,x) defined on the spacetime manifold of the Einstein universe. 

Now suppose that there is an abstract Hilbert space HI of vectors |/). Let x and t be the operators of the space 
and time coordinates, correspondingly, defined on H. Moreover, let \t, x) = \x) be the complete set, or basis, of 
eigenvectors of x and t, i.e. x|x) = x\x), t\x) = t\x). The set |x) is usualy called the coordinate basis in EL Obviously, 
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the completeness and normalization conditions for the coordinate basis |x) are valid: 

d 4 x % /=g\x)(x\=I, (Al) 

(x'\x) = 5{ ~ x ~^\ (A2) 

where I is the unit operator in EI and g = det g^ v . Due to (|A1[) , it is possible to expand any vector |/) £ H in terms of 
the basis |x), namely |/) = J d 4 x^/—g\x)(x\f). The quantity (x\f) is called x- (or coordinate) representation for the 
vector |/). Identifying (x\f) with f(t,x), we see that the Hilbert space H of all quadratically integrable functions is 
simply the coordinate representation of the above introduced abstract Hilbert space EL Furthermore, if in the Hilbert 
space H there is an arbitrary operator a, then the matrix A, whose matrix elements in the coordinate basis are just the 
quantities (x'|a|a;), is called the x- (or coordinate) representation of the operator a. Obviously in order to define an 
operator a in the abstract Hilbert space H, it is sufficient to define its matrix A = (x'\a\x), which acts in the Hilbert 
space H of all quadratically integrable functions, in the coordinate basis. It is clear that TiA =J d 4 x^/—g(x\a\x). 

Now, let us consider in H the two commuting operators p»o and h such that in the ^-representation they look like 
Po = (x'\po\x) and H = (x'\h\x), correspondingly (see section HTT1) . The operators Pq and 7i have a common set of 
eigenfunctions ^i ar ip (t, x) defined in section [TTT1 from which it follows that 

d i Xy/^g^ lm)po (t,x)% a , v , p , o (t,x) = 2TrS(p' a - p )Sii'S aa >S, ir) > . (A3) 

The eigenfunctions ^i arj p (t, x) are the coordinates of the corresponding eigenvectors \larjpo) £ H of the operators po 
and h (recall, / = 0, ...oo; a — I, .., di = 2(1 + l)(l + 2); -q = ±1; — oo < po < oo): 

\larjpo) = / d 4 xy/^g\x) 1 $i a , 1Po (t, x) = j d i x^/ zr g\x)(x\la-qp a ). (A4) 

(Clearly, po|/ar/po) = Po\l a VPo} an d h\lai]po) — r]Ei\larjpo}.) It follows from (|A4|) that ^i a rip (t, x) = (x\lar)po)- 
Using this relation in the normalization condition (|A3[) and then integrating there over x with the help of (|Aip . we 
obtain 

(l'a'r)'p' \laripo) = 2ir6(p' - p Q )6u>5 aa >6 rn > (A5) 



which is the analogue of the normalization condition (|A2|) . It is possible to show that the completness condition for 
the basis \laripo) follows from (|A5[) : 

^2 / -^-\ lar JPo)(la'npo\ = I. (A6) 

lar/ 

Let us construct in EI the following operator 

6 = A 2 - (p + m) 2 + (h + Stf (A7) 

which is diagonal in the basis (|A4[) . i.e. each vector \larjpo) is its eigenvector with corresponding eigenvalue £i avPo 
In the coordinate representation its matrix (a/|6|x) coincides with the operator O PD)) . So, 

Tr6= f d 4 xV^{x\6\x) = Jd^x^-gY, J ^ E / ^{x\lar,p ){lar 1 p \6\l'a'r ] 'p l )(l'a'r, l p l o\x), (A8) 

lar] V a.' rj' 



where the last equality was obtained by employing the completeness relation (|A6[) . Now, by using in this formula the 
eigenvalue condition d\lanp ) — Siocnpol^VPo) j t ne normalization condition (|A5|) . and, finally, by performing in the 
obtained expression the integration and summation over primed indices, we have 

Tr6= / d 4 xv r ^g~^2 / -^-£i a vpo( x \ la VPo)(lar]p \x} = / d 4 xy/=g~^2 / l^^rip a ^i am y a {t,x)^* lanpo (t,x). (A9) 

lar] ' lar] 

Since in (|A9[) the quantities yf—g and 1 $i ariP0 {t, x)^* (t, x) =ipi ari (x)ipiari{x) (the last relation is due to the notations 
accepted in formula (|2"2"|) and below) do not depend on the time coordinate, the expression (|A9|) is proportional to 
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the infinite time interval T = J dt. The remaining x-integration in (|A9|) gives simply unity due to the relation (1221) . 
So we have 

TrO = T^ j ^ - T E / |f * ^ " (P0 + ^ + ^ + ^ ' (A10) 

where the fact that each eigenvalue Si ar)Po is (i/-fold degenerated is taken into account and the notations from (|21l) - (|2"5)l 
are used. In a similar way it is possible to obtain the quantity TrlnO: 

Tr In d = Tj2J^ k ^Wo = T E / f^ ln t A ' " &* + ^ + ^ + ^] • ( An ) 

APPENDIX B: THE ROLE OF QUANTUM FLUCTUATIONS AND FINITE SIZE EFFECTS 

It is well known that spontaneous symmetry breaking in low dimensional quantum field theories may become 
impossible due to strong quantum fluctuations of fields. The same is also true for systems that occupy a limited 
space volume. However, as it is clear from physical considerations, the finite size in itself may not practically forbid 
the spontaneous symmetry breaking, if the characteristic length of the region of space occupied by the system is 
much greater than the Compton wavelength of the excitation responsible for tunneling and restoration of symmetry. 
(Indeed, one may recall here well known physical phenomena such as the superfluidity of Helium or superconductivity 
of metals that are observed in samples of finite volume) . This idea has been discussed for some scalar field theories in 
the Einstein universe for instance, in [40, |4l[. In this Appendix, we shall demonstrate that, under certain conditions, 
dynamical symmetry breaking in NJL-type models is indeed possible in the closed Einstein universe. In particular, 
we will show that, if the radius of the universe is large enough such that the fluctuations of quantum fields are 
comparatively small, the symmetry breaking obtained in the mean field approximation is not forbidden. 

For illustrations, let us confine to the analogous case of the chiral condensate and consider the simplified case of 
the linearized Lagrangian © with /x = 0, S/j, = 0, m = 0, 7Tfc = 0, 

I = q{i^ v V v - a)q - ^^ 2 » 
and the corresponding partition function Z = e ,ff with 

/2 
d*xV=g?g - ?Trhi(yfV v - <r). (Bl) 

Now, supposing that a = ctq + <f>, where <7o is the vacuum expectation value of the field a and <f> denotes its quantum 
fluctuation, we obtain: 

ln(P - <f>) w ln£> - £> _1 (/> - -{V 1 0)(X> _1 </>) 

where V = vf^J v — uq. Thus Z = ZqZ^, where Zq = exp iSo, So = — J d^x^J—gVo, and 

y ° = 4t + 7^v=f rlnP (B2) 

is the effective potential at a = ero- The contribution of quantum fluctuations up to the </> 2 -term to the effective action 
S^ is given by 

i>exp(iS^) = Jtytsq? l-i j dtxyTg-^tt 2 + W) - Tr(p-^ + ^(P" 1 ^" 1 *)) 

It is evident, that in the above expansion the term linear in (f> corresponds to the so-called tadpole diagram with one 
external 0-line and the term quadratic in <f> corresponds to the "polarization operator" diagram of the <j> field with 
one fermion loop. From (|B2[) we can write the stationarity condition and find the gap equation, dV^/daa = 0, 

^jd i x^—g = iTvV-\ (B3) 



and the linear terms in 0, corresponding to the tadpole diagram, cancel out in Z^. 
fluctuations to the field action is given by 
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Thus, the contribution of 



Sd 



d^V~9^ + ^(V-^)(V-\ 



(B4) 



Next, we shall calculate the contribution of fluctuations cf> to the effective action, taking into account the quark loop in 
the gravitational field. (Note that this corresponds to the integral (11) and Fig. la in [lj].) For our purpose of making 
estimates of the role of fluctuations, it is sufficient to limit ourselves to the consideration of fluctuations depending 
only on time. In this case, we can extract the necessary kinematical factor for the meson fluctuation field and then 
integrate over the time-component of the loop momentum in the limit of vanishing external momenta. Finally, after 
going to the basis for the Dirac equation in the Einstein universe (see (J2T))) - (|2"2"T) ) we obtain a sum over fermion loop 
quantum numbers I instead of an integration over momenta of free quarks made in [l| . The sum is divergent and we 
regularize it by the cut off A. In this way we obtain the effective action 



S,i, 



d 4 x. 



-9 



£_ 
AG 



— — y, 

47r 2 a 3 Z-~< ' 

1=0 



</ A 2(Z + l)(Z + 2) 



Ei 



Aafa 2 



{dtcpf 



AEf 



(B5) 



where Ei and uii are given in (|2ip . The summation over I in the first term in parenthesis in the above equation cancels 
out by the term cj) 2 /AG, due to the stationarity condition (|B3j) 



1 = 



3G 



1=0 



111 + W + 2)±. 



(B6) 



After this we obtain 



Sa — 



3 1/1 

dt- X; e~^ A 2{l + l)(l + 2)^ -tf olc 

i=o l ^ 



or S,/, — J dtLcf, with the Lagrange function 



!(• 



4cr^ 2 ] VZ- 1 



where V — j d x^/—g = 2n 2 a 3 is the space volume, and the renormalization 2-factor is defined as 



4V<^ 



e -"i/*2(l + l)(l 



K 2 + ^o 2 ) 3/2 ' 



(B7) 



(B8) 



For comparision, let us consider the limiting case of flat space which is reached by the replacements 

dl, (l + l)(l + 2)^fa 2 



E 
i 



in (|B8[) . Then the Z- factor takes the form of the integral in the Euclidean spacetime 

d A p 1 
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{2-Kfipl + eir 



with £ p being the quark energy. This expression evidently corresponds to the similar formula for the Z-factor in the flat 
space case considered in (lj. Next, let us perform the field renormalization (j> = Z^^> r in (|B7|) . Quantum fluctuations 
of the boson field near the ground state ag can now be estimated, if we consider the renormalized expression (|B7|) as 
the Lagrange function for a harmonic oscillator (here, we follow the idea of [40] 3 ) with the mass m and frequency w, 
formally given here by the relations 



to = V, uj 2 = Aa 2 = M 2 . 



3 In our case of the NJL model, the consideration of the quark loop diagram is essential (see [lj). This differs from Ref. |4d |. where the 
</> 4 — model of a self-interacting scalar field was considered and the scalar loop contribution to the fluctuation Lagrangian was calculated. 
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Then we can estimate the quantum fluctuations as 

1 1 1 



' ' raw V2ct VM 
where M is the mass of the composite er-meson. Thus, we obtain 



P)*«< (B9) 

The estimate (|B9|) gives a criterion for the role of quantum fluctuations for a system with finite volume. Clearly, 
quantum fluctuations can be considered negligible, if Oq 3> (4> 2 )- This is in agreement with the physical requirement 
that quantum fluctuations should be negligible if aao 3> 1, i.e., if the radius of the universe is much greater than the 
Compton wavelength A = jj = ^- of the cr-meson (quarks) (see also [40| ) 4 . 

The above estimates are certainly of a qualitative nature, and hence they do not allow us to find an exact value 
of a critical radius such that symmetry breaking for lower values of the curvature radius of the Einstein universe is 
forbidden. 
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